Liquid-like behavior of supercritical fluids 
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The high frequency dynamics of fluid oxygen have been investigated by Inelastic X-ray Scattering. 
In spite of the markedly supercritical conditions (T ~ 2T C , P > 10 2 P C ), the sound velocity exceeds 
the hydrodynamic value of about 20%, a feature which is the fingerprint of liquid-like dynamics. The 
comparison of the present results with literature data obtained in several fluids allow us to identify 
the extrapolation of the liquid vapor-coexistence line in the (P/P c , T/T c ) plane as the relevant 
edge between liquid- and gas-like dynamics. More interestingly, this extrapolation is very close to 
the non metal-metal transition in hot dense fluids, at pressure and temperature values as obtained 
by shock wave experiments. This result points to the existence of a connection between structural 
modifications and transport properties in dense fluids. 

PACS numbers: 67.55. Jd,67.40.Fd,61.10.Eq,62.50.+p 



With the recent advent of the Inelastic X-Ray tech- 
nique (IXS) [l|, 0] , that complemented the Inelastic Neu- 
tron Scattering one extending the accessible exchanged 
momentum Q to the whole first pseudo Brillouin zone, 
the collective dynamics of disordered systems have been 
investigated in different classes of liquids and glasses 
[!, 0, Among other features, the specific interest 

here is the discovery that in all the investigated liquids 
and glasses it exists a positive dispersion of the sound 
speed, i.e. an increase of the velocity of the longitudi- 
nal acoustic mode from the (relaxed) ultrasonic value to 
an infinite frequency (unrelaxed) value. Such anomalous 
dispersion is ascribed to the presence of one (or more) re- 
laxation processes interacting with the dynamics of the 
density fluctuations. 

On a qualitative ground, a relaxation process is char- 
acterized by a specific time, r, which marks the border 
between "viscous liquid" -like dynamics (for lot « I) to 
"elastic solid" -like one (for cut >> I). Correspondingly, 
the longitudinal sound velocity, i. e. the velocity of prop- 
agation of the density fluctuation, undergoes a transition 
(positive dispersion) from its "low" frequency limit, Co, 
which characterizes the "liquid" value, to its "infinite" 
frequency limit Coo > Co, characteristic of the "solid" re- 
sponse of the system. 

Usually, more than one single relaxation process is ac- 
tive, and a " sequence" of sound velocity dispersions takes 
place. In particular, convincing evidence has been re- 
ported in the last years supporting a scenario in which 
the density fluctuations decay according to three differ- 
ent relaxation mechanisms [g]: i) a thermal process (t^), 
arising from the coupling of temperature and density, 
which marks the transition from adiabatic (iOT t h » 1) 
to isothermal (ivTth << 1) dynamics, ii) a structural 
(or a) process, strongly temperature dependent, whose 
timescale r Q is directly related to the viscosity and hence 
to the mass diffusion coefficient, in) a microscopic pro- 



cess, which drives the decay of the density fluctuations 
at a given wavelength, induced by the underlying disor- 
dered structure that does not support this wavelength 
as an energy eigenstate Q. This scheme, validated by 
several experimental and numerical evidences, has been 
developed for glasses and normal liquids, i.e. for systems 
well below their liquid-gas critical point. 

In this letter we investigate whether this approach can 
be applied to the supercritical phase, i.e. if and how the 
dynamics of a dense liquid evolves on crossing the crit- 
ical temperature. On a qualitative ground, one might 
expect that on abandoning the liquid phase the posi- 
tive dispersion should vanish. As matter of fact both 
the structural and microscopic processes are likely to be- 
come inactive as soon as the dynamics loses its coopera- 
tive nature and the concept of instantaneous vibrations 
is no longer well defined. Looking at the available IXS 
literature data for simple liquids and supercritical fluids 
d EE EE EI ES El EI ES El, one may actually con- 
clude that the positive dispersion vanishes on approach- 
ing the critical temperature. However, this simplified 
view seems to be contradicted by the results reported 
here for supercritical Oxygen: at twice the critical tem- 
perature we still find a 20 % positive dispersion. This 
result seems to indicate that the end of the liquid dy- 
namics is not only marked by the critical temperature. 
Moreover, by examining the whole set of existing data, 
we conclude that the extension of the liquid-vapor coex- 
istence line into the supercritical region (corresponding 
to the Widom's line, i.e. the thermodynamic locus of the 
specific heat maxima) defines the transition from liquid- 
like to gas-like dynamics. This finding could help to give 
insights to the phenomena of the non metal- metal transi- 
tion in hot dense simple fluids such as Hydrogen, Oxygen 
and Nitrogen, that has been observed at thermodynamic 
points which actually lie on the extrapolation of the co- 
existence line. 
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The experiment was carried out at the IXS Beamline II 
(ID28) at the ESRF. The X-rays from an undulator are 
monochromatized by a double crystal monochromator 
and a high-energy resolution backscattering monochro- 
mator, operating at the silicon (9 9 9) reflection order. 
The backscattered photons of energy 17794 eV (3xl0 10 
photon/s) are focused by a cylindrical mirror and a 
Kirkpatrick-Baez multiplayer, providing a focal spot of 
30 (horizontal) by 80 (vertical) ^m 2 Full Width at Half 
Maximum (FWHM) at the sample position. The scat- 
tered photons are energy-analyzed by a Rowland circle 
five-crystal spectrometer, with an overall energy resolu- 
tion 5E « 3.0 meV, and detected by a Peltier-cooled sil- 
icon diode. The momentum transfer Q = 2ki ■ sin(0 s /2), 
where ki is the incident photon wave vector and S is the 
scattering angle, is selected by rotating the spectrome- 
ter around a vertical axis passing through the scattering 
sample. Spectra at five different momentum transfers 
can be recorded simultaneously by five independent an- 
alyzer systems with a resolution SQ « 0.2 nm _1 . We 
have used a diamond anvil cell (DAC) of the membrane 
type, equipped with 600 /im culet diamonds. The cell 
was loaded by condensing high purity liquid oxygen in a 
sealed vessel at low temperature. Pressure is measured 
by the shift of the ruby fluorescence wavelength [3l| . The 
sample dimensions were about 300 /im in diameter, and 
80 /xm in thickness. The IXS measurements were per- 
formed at room temperature at 0.88, 2.88 and 5.35 GPa 
(±0.01 GPa). The highest pressure experimental point 
is close to the melting line, but still in the fluid phase as 
was obtained by melting the crystal along an isothermal 
decompression (300 K) under visual inspection by using 
a microscope. 

The evolution of the IXS spectra as a function of the 
momentum transfer, Q, is shown in Fig[T]for the 3.0-^12.6 
imr 1 Q-range and for the three investigated pressure 
values. The inelastic signal can be clearly observed at 
the two sides of the quasi-elastic line: it is structured at 
the low Q-values, and it becomes broad and featureless 
on increasing Q, although it is still well visible on top of 
the tails of the quasi-elastic peak. 

The IXS spectra are related to the classical dynamic 
structure factor S(Q,u) through the detailed balance 
factor and the instrumental resolution function, R(lu), 
through: 
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FIG. 1: Selected IXS spectra of supercritical oxygen at T=300 
K. The three columns correspond to three different pressure 
values (0.88, 2.88 and 5.35 GPa from left to right) and the 
different rows report spectra taken at the indicated Q values 
(3.0, 5.4, 7.8, 10.2 and 12.6 nm" 1 from bottom to top). 
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In this expression ujq = kgTQ- /mS(Q), being m the 
molecular mass, kg the Boltzmann constant and S(Q) 
the static structure factor, and M(Q,ui) = M'(Q,u>) + 
iM"(Q,uj) is the memory function in the frequency do- 
main. In the hydrodynamic limit, the decay of the mem- 
ory function proceeds via a timescale related to the ther- 
mal diffusivity Tth — "fDxQ 2 , where 7 is the specific heat 
ratio. A second decay, whose timescale (r a ) is related 
to the typical structural relaxation time, is introduced 
to account for the viscoelastic nature of the system. Fi- 
nally, a third, microscopic, relaxation time is introduced 
to represent the effect of the disordered instantaneous 
atomic arrangement on the propagation of the density 
fluctuation at a given Q. As in most cases this third 
contribution can be considered as " instantaneous" , and 
represented by a delta- function. Consequently, the mem- 
ory function in the time domain reads: 



A model for the classical S{Q,iv) has been obtained M(Q,t) = (7 - l)w 2 e~ 7l>rQ2t + 2TQ 2 5{t) + A 2 L e' t/T " 

exploiting the generalized hydrodynamic approach, in (3) 

which the dynamic structure factor is related to the sec- In this expression TQ 2 is the amplitude of the micro- 

ond order memory function through (l8| : scopic process, while A 2 = (c 2 ^ — c 2 )Q 2 is the intensity 
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of the structural process, being c and Coo the low and 
high frequency sound velocities respectively . By plug- 
ging this expression into Eq. @, and using Eq. ([TJ) , we 
were able to describe the experimental data, as shown by 
the best fit reported in Fig. [TJ In the fitting procedure, 
the parameters A, r a , T and luq were left free, Dt and 7 
were fixed to their Q=0 values obtained by the Impulsive 

3 
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FIG. 2: Dispersion relation of supercritical oxygen at T=300 
K at the three measured pressures (circles: P=5.35 GPa, tri- 
angles: P=2.88 GPa, diamonds: P=0.88 GPa), from the cur- 
rent maxima associated to the best fit according to Eq. J2). 
The dotted lines indicate the adiabatic sound velocities as 
determined by low frequency (GHz) ISTS measurements 
while the continuous lines are a sine- like best fit whose Q — > 
slope represents the high frequency sound velocity. 



P [GPa]) cisTs\m/s\ c IX s\m/s] 
0.88 1920 2340 

2.88 2980 3600 

5.35 3680 4440 



TABLE I: Low and high frequency sound velocity in super- 
critical oxygen. Experimental uncertainties are ±50 m/s. 

In Fig. [2] we report the dispersion curves as deter- 
mined by IXS at the three investigated pressures. These 
were obtained from the maxima of the current correlation 
spectra J(Q,w) = co 2 /Q 2 S(Q, lj) (circles), utilizing the 
fit results for the model S(Q,u>). For space reason the 
values of the other fitting parameters are not reported 
here. For each pressure, the low Q slope of the evolu- 
tion of the energy values represents the sound velocity 
measured at high frequency. We note that the observed 
high frequency sound velocity values always lie above 
the low frequency, adiabatic sound velocity measured by 



the ISTS technique [19|(see Table 1) [32j. This result 
demonstrates the presence of positive dispersion -the fin- 
gerprint of a liquid-like dynamics- in a supercritical fluid 
well above its T c and P c . As an immediate consequence, 
one has to rule out the idea that the critical temperature 
alone marks the boundary of simple liquid dynamics. 

In order to asses what really determines the end of the 
simple liquid dynamics regime, we compared the present 
results for oxygen to existing literature data. In fig [3] we 
highlight the presence of positive dispersion in a library 
of fluids at several thermodynamic points. Full and open 
circles indicate a positive dispersion above and below 5%, 
respectively. As it can be seen, systems around or mod- 
erately below the critical temperature, but well above P c 
still exhibit positive dispersion, which fades out on in- 
creasing temperature. What is really surprising is that 
positive dispersion still persists in a strongly supercriti- 
cal system, the case of Oxygen reported here, in which 
both temperature and pressure are well above the critical 
point. In order to rationalize such a puzzling behaviour, 
we report in F ig. JH the Plank- Riedel equation for the co- 



existence line [21|: 
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The values of a, b and c (a=4.270, b=-4.271 and 
c=1.414) have been obtained by fitting the liquid-vapour 
coexistence lines of Neon, Oxygen and Nitrogen [22j ]. 
Eq. |4] can represent the noble gases liquid-vapour co- 
existence line 22j in the spirit of the principle of cor- 
responding states, but it also describes rather well the 
liquid-vapour coexistence line for the other simple fluids 
discussed here (see Fig. [3]). This curve is extrapolated 
here above the critical point. On the basis of Fig. [3l one 
is led to conclude that the end of simple liquid dynamics 
is identified by the extrapolation of the coexistence line 
although, on a strict thermodynamical ground, such line 
has a well defined meaning only below the critical point. 

Recent shock wave measurements have shown a density 
and temperature driven non metal-metal transition in 
hot dense fluid Oxygen [23| , Nitrogen Q and Hydrogen 
[25| at thermodynamic points which surprisingly lie on 
the proposed extrapolation line, much above the region 
where acoustic properties have been investigated (inset of 
figure [3]). Also in the case of Mercury, which is a metallic 
liquid and an insulator gas 26], the non metal- metal 
transition in the supercritical fluid, observed slightly 
above the critical point, takes place around the extrap- 
olation of the liquid-vapour coexistence line, see Fig. [3] 
As a matter of fact, the metal-non metal transition of ex- 
panded mercury, which has been extensively studied 27 1 
and references therein) , appears to be related to the den- 
sity decrease which takes place on crossing the extrapo- 
lation of the liquid-vapour coexistence line. Furthermore 
the positive dispersion disappears on crossing the same 
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extrapolation, as shown in Fig. [3l so that in this case we 
have the exact coincidence of the two phenomena: the 
change in dynamics and the metallization transition. 
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simple fluids into a gas-like and a liquid-like domain, re- 
spectively, which differ in the local structure, resembling 
the sub-critical behaviors. Specifically, the "liquid-like 
phase" is expected to be denser and less rich in entropy 
than the "gas- like phase" . These differences, even if sub- 
tle at a first glance, could support the major changes 
observed in the dynamical and transport properties. If 
this point will be confirmed by further experiments on 
different fluids and on an extended thermodynamic re- 
gion, it would imply a profound interconnection between 
structural modifications in fluid materials occurring at 
very different P-T conditions which, in turn, could stim- 
ulate and support a new extended and unified view of 
fluid matter thermodynamics. 

We acknowledge Alexandre Beraud for assisting us 
during the experiment. 



FIG. 3: Sketch of the (P/P c -T/T c ) plane. The black lines 
are the best fit of the average of the liquid-vapour coexis- 
tence lines for Neon, Nitrogen and Oxygen using Eq. [4] (see 
text), drawn below (continuous line) and extrapolated above 
(dashed line) the critical point, respectively. The dashed blue 
line is the melting line of Oxygen [22l |28{ . Each color corre- 
sponds to a different investigated system. Open points repre- 
sent cases where the positive dispersion of the sound velocity 
has not been observed, full points cases where there is a clear 
signature of positive dispersion. Data for liquid metals are 
taken from [3] and references therein (see Table I); data on 
Neon are from [2S| and from [l5l . Il6l ]; data on Mercury are 
from [27l ]; data on Nitrogen are from [l7j ] . Inset: thermody- 
namic points where the non metal-metal transition has been 
experimentally observed for Mercury (ref. [26|), Oxygen (ref. 
[13|), Nitrogen (ref. [Hj]), and Hydrogen (ref. The last 

three cases were investigated by means of shock wave experi- 
ments. 



In conclusion, we have shown here that a strongly su- 
percritical fluid (Oxygen at T/T c w 2 and P/P c > 100) 
exhibits dynamical properties that are commonly as- 
cribed to ordinary liquids. Based on the comparison with 
the behavior of several other fluids as a function of their 
thermodynamic state, we conclude that the extrapola- 
tion well above the critical point of the liquid-vapor co- 
existence line marks the boundary between simple liquid 
dynamics and the collision dominated regime character- 
izing the gas phase. Furthermore, the existing data on 
the location of the non metal- metal transition of mercury 
and of hot dense simple fluids obtained by shock wave 
experiments indicate that the structural changes which 
take place on crossing the extrapolation of the liquid- 
vapor coexistence line are also responsible for the non 
metal-metal transition in overcritical fluids. 

We suggest that the extension of the coexistence line 
beyond the critical point splits the P-T phase diagram of 
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